Binary Nonlinearization of Lax pairs of Kaup-Newell Soliton Hierarchy 

Wen-Xiu Ma^, Qing Ding"!", Wei-Guo Zhang* and Bao-Qun Lut 

^Institute of Mathematics, Fudan University, Shanghai 200433, P. R. of China 

^FB Mathematik und Informatik, Universitdt Paderborn, D-33098 Paderborn, Germany 

*Dept. of Math. Sz Mech., Changsha Railway University, Changsha 410075, P. R. of China 



Summary. — Kaup-Newell soliton hierarchy is derived from a kind of Lax pairs 
\Q ' different from the original ones. Binary nonlinearization procedure corresponding to the 

| Bargmann symmetry constraint is carried out for those Lax pairs. The proposed Lax 

pairs together with adjoint Lax pairs are constrained as a hierarchy of commutative, finite 
dimensional integrable Hamiltonian systems in the Liouville sense, which also provides 
us with new examples of finite dimensional integrable Hamiltonian systems. A sort of 
involutive solutions to the Kaup-Newell hierarchy are exhibited through the obtained 
finite dimensional integrable systems and the general involutive system engendered by 
^T) . binary nonlinearization is reduced to a specific involutive system generated by mono- 

o 
o 

00 . 

, PACS 03.40K - waves and wave propagation: general mathematical aspects 

1. Introduction 

It is very important (and interesting) to search for integrable reductions of soliton equations. 
Similarity reductions lead to equations of Painleve type [1], and algebraic pole expansion solutions 
to soliton equations engender finite dimensional integrable A^-body Hamiltonian systems [2]. These 



X 



nonlinearization. 



are two main examples of integrable reductions. 

Recently a systematic approach called the nonlinearization (more precisely, we call it mono- 
nonlinearization) of Lax pairs is introduced for generating finite dimensional completely integrable 
systems from soliton hierarchies by Cao and Geng et al. [3-5]. Based upon the idea of nonlin- 
earization, one of the authors proposed binary nonlinearization technique for soliton hierarchies 
[6,7], which was successfully applied to a few well-known soliton hierarchies [6,7]. Nonlinearization 
method may yield a large variety of finite dimensional completely integrable Hamiltonian systems 
whose independent integrals of motion can be explicitly given out, and it provides a way of solving 
soliton equations by separation of spatial and temporal variables. 

The integrable reductions in the mono-nonlinearization and binary nonlinearization are related 
to a kind of specific symmetry constraints [6,7]. Bargmann symmetry constraints (usually the first 
order with respect to derivative) lead to the most interesting reductions in the nonlinearization 
because of their simplicity. These constraints require conserved covariants to be linear functions 
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but not differential functions with respect to the potentials. The higher order symmetry constraints 
need to introduce new dependent variables, the so-called Jacobi-Ostrogradsky coordinates. The 
spatial parts of the resulting nonlinearized Lax pairs by nonlinearization are exactly the restricted 
flows named by Antonowicz and Wojciechowski etc. [8]. 

In this paper, we would like to discuss binary nonlinearization for Kaup-Newell hierarchy [9] of 
soliton equations. In Section 2, following the standard procedure, we shall reconstruct Kaup-Newell 
soliton hierarchy from a particular spectral problem other than the original Kaup-Newell spectral 
problem [9] and give rise to some algebraic properties connected with Kaup-Newell hierarchy. 
In Section 3, we shall in detail exhibit binary nonlinearization procedure corresponding to the 
Bargmann constraint with zero boundary conditions and without zero boundary conditions and 
vigorously constrain the proposed Lax pairs and adjoint Lax pairs of Kaup-Newell hierarchy into 
a hierarchy of commutative, finite dimensional integrable Hamiltonian systems in the Liouville 
sense. By the resulting finite dimensional integrable system, we shall further exhibit a sort of 
involutive solutions to Kaup-Newell hierarchy. Finally in Section 4, we reduce the general involutive 
system derived from binary nonlinearization to a specific involutive system generated by mono- 
nonlinearization and give some other remarks. 

2. New Lax pairs of Kaup-Newell hierarchy and some related properties 

Let us consider the following specific spectral problem 

<f> x = U<f> = U(u,X)<f> = ( t \)<f>, <t>=(t 1 ), u=(l). (2.1) 



^ Ar —A I ' V (j>2 J ' \r 
This spectral problem is a little different from the original Kaup-Newell spectral problem [9]. We 
find by a large amount of computation that it is a possible spectral problem leading to an integrable 
hierarchy. Although we don't know if there exists a gauge transform between (2.1) and the original 
spectral problem, we shall see later that its related integrable systems may be transformed into 
Kaup-Newell hierarchy. 

In order to derive isospectral systems associated with (2.1), similar to Ref. [10] we first solve 

the adjoint representation equation (see Ref. [11]) V x = [U,V] of <p x = U<p. Take 

a A"V 
c —a 



V 

Noting that 



[U,V] 



( qc — rb 2b — 2qa 
\ 2\ra — 2Ac rb — qc 



we find that the adjoint representation equation V x = [U, V] becomes 

' a x = qc- rb, 



X~ 1 b x = 2b- 2qa, (2.2) 
c x = 2\ra — 2Ac. 
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On setting a = £\> ctjA l , 6 = X^i>o M *> c = X^>o ( 2 - 2 ) y ields equivalently 

' 26 — 2gao = 2rao — 2co = 0, 

tt^x 

b ix = 2b i+ i - 2qa i+1 , i > 0. 
, c ix = 2ra i+1 - 2c i+ i, 



(2.3) 



We choose ao = 1, bo = q, cq = r and assume that ai|[ u ] = o = &i|[ u ]=o = c i|[u]=o = 0, i > 1) where 
[u] = (u, u x , ■ ■ ■ ) (this means to select zero constants for integration). Because we have 

2a i+ i iX = 2qc i+1 - 2rb i+1 
= q(2ra i+1 - c ix ) - r(b ix + 2qa i+1 ) 
QCix fbi x , 

the equality (2.3) gives rise to a recursion relation for determining a,, bi, q: 



a i+ i = ~-d~ l (qc ix +rb ix ), 

b i+i = 2 bix ~ 2 qd l ( qCix + rbix) ' i ~ °' 
Ci+i = -^ c in; - ]^rd~ 1 {qc ix + rb ix ), 



Therefore, for instance, we can obtain 



(2.4) 



fl2 = ^q 2 r 2 + ^gr x - ^q x r, b 2 = ^q xx - ^qq x r + ^q 3 r 2 , c 2 = ^ „ ; 



From (V 2 ^ = [£/,y 2 ], we see that (\trV 2 ) x = (a 2 + A" 1 ^ = 0. Thus we have a 2 + A _1 6c = 1 
by observing (^trV 2 )|[ u ] =0 = 1. Further we may obtain 



m — 1 



m— 1 



2a r 



X! a » a ™-' ~ X] fe i C "i-l-i 5 "I > 2. 



(2.5) 



?=1 



i=0 



This equality shows by an argument of mathematical induction that en , 6j , q are all differential 
polynomial functions of u. 

Now we introduce the following auxiliary problem associated with the spectral problem (2.1) 



<p tm = v^U, = (\ m+1 v) + - ( 



C m + 1 — Om+1 



A(A m a)+ (X m b) + \ 
A(A m c)+ -A(A m a) +y l ' m - U ' 

(2.6) 
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where the symbol + stands for the selection of the polynomial part of A. Then cross differentiation 
(i.e. taking d/dt m (2.1), d/dx (2.6) and setting them equal) leads to isospectral integrable systems 

«*» = (;) = Km = (In) = * mux > m -°- (2 - 7) 

Here the operator $ reads 

\ -\drd~ x r -\d-\drd~ x q)' °* 
which may be verified to be a hereditary symmetry operator [12], i.e. to satisfy 

<f> 2 [K, S] + [&K, $S] - §{[K, $S] + [&K, S}} = 0, 

for arbitrary vector fields K,S. Obviously, under ■£ — > — iJ^, ^| > —i^—, (2.7) is transformed 

into Kaup-Newell hierarchy [9]. We call (2.7) as a kind of real form of Kaup-Newell hierarchy 
because (2.7) itself is a hierarchy of real systems. The first two nonlinear systems in the hierarchy 
(2.7) are as follows 
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9ti = ^ [Qxx - (q 2 r) x ] , 



r tl = -\[rxx + (qr 2 ) x ]; 



n 2 = ~£xxx + j(qrr x ) x + ^{q 2 r 3 ) 



2 

2 

The former system may be rewritten as the derivative nonlinear Schrbdinger equation [9] 

1 1 

nu = -^q xx + ^i(q*q 2 ) x 

under £ -> ^ -> -i^-, 9 = -r*. 

We below want to tersely exhibit bi-Hamiltonian structures of Kaup-Newell hierarchy (2.7) by 
applying a powerful tool, i.e. the so-called trace identity proposed in Ref. [13,14]. Towards this 
end, we need, as is usual, the following quantities which are easy to work out: 

<V,— >=2a + \ 1 rb, < V, — — >= c, < V, — — >= 6, 
aA or 

where < •, • > denotes the Killing form: < A, B >= tr(AB). Then by means of the trace identity 
[13] 

5 Tr dU ^ d Tr dU Tr dU , T 
5u < 9X >= A_1 9A A1(< * >• < V > 8F » ' ■» = c ° nst " 

we obtain immediately 

A(2a + A- 1 r6)=A-^A-(c,6) T . 



Comparing the coefficients of A m 1 on two sides of the above equality and noticing the equality 
(2.3), we have 

(5 1 & 

— [2a m+1 + -{rb m + qc m )] = — (2a m+1 + rb m ) = (7 - m)(c m , b m ) T , m > 0. (2.8) 

The equality (2.8) with m = yields that the constant 7 = 0. Therefore (2.8) gives rise to an 
important formula 

5 1 

— [2a m+ i + -{rb m + <?c m )] = -m(c m , b m ) T ', m > 0. (2.9) 
At this stage, Kaup-Newell hierarchy (2.7) may be transformed into the Hamiltonian form 

"'- = ('),. = A ™ = ( tl) - J0 ™ - J *™ (*) - m 2 °- (2 - 10> 

Here the Hamiltonian operator J, the recursive operator ^ and the Hamiltonian functions H m are 
defined by 











: 


, * = $* = ( 









-\d — \rd 1 qd —\fd l rd 
— \qd~ 1 qd \d — \qd~ 1 rd ^ 



(2.11a) 



H = -2a!, H m = — -\2a m +i + \{rb m + qc m )], m > 1. (2.11b) 
m L 2 J 

A direct calculation can show that J and J*S> constitutes a Hamiltonian pair. Therefore the flows 
determined by Kaup-Newell hierarchy commute with each other, and every Kaup-Newell system 
has an infinite number of symmetries {K n }^ =0 and conserved quantities {H n }^ =0 . It is known that 
Magri [15] showed the existence of bi-Hamiltonian formulation of Kaup-Newell hierarchy, however 
here we also give the explicit expressions (2.11b) of the Hamiltonian functions H m . 
Further, similar to Ref. [16], we can get the Zakharov-Shabat equation 

Vt ] ~ Vtt ] + [V {m) ,V^] = (V (m ^'[K n ] - (V^)'[K m ] + [V^ m \V^] = 0, 

which also shows the commutability of the flows of (2.10). In addition, we can directly verify 

V tm = [V {m \Vl m>0, (2.12) 

whenu tm =K m ,i.e. U tm -V^ m) + [U, V^] =0, m > 0. In fact, we may deduce that V tm -[V^ m \ V] 
satisfies the adjoint representation equation of <p x = U<f> and that (V tm — [y( m ), V])|[ u ] = o = 0. Hence 
(2.12) follows from the uniqueness property [6] (i.e. if V x = [U,V] and V|[ u ] = o = 0, then V itself 
vanishes) of the adjoint representation equation. The equality (2.12) shall be used to elucidate the 
commutability of the flows generated by nonlinearization in the next section. 
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3. Binary nonlinearization of Lax pairs 

In this section, we would like to exhibit a binary nonlinearization procedure, as in Refs. [6,7], 
for the Lax pairs and adjoint Lax pairs of Kaup-Newell hierarchy (2.7) or (2.10). Note that U tm — 
+ V^] = if and only if {-U T ) tm - {-{V^) T ) X + [-U T ', ~{V^) T ] = 0. Kaup-Newell 
systems (2.10) have another kind of Lax pairs called the adjoint Lax pairs 

U x = -U T ip = -U T (u, \)il>, (3.1a) 
\ A m = -{V {m) ) T 4> = -(V^) T (u, X)iP, (3.16) 
where T means the transpose of matrix and ip = (^1,^2) ■ The Lax pairs and adjoint Lax pairs 
are a basic object of binary nonlinearization. 

3.1. Bargmann symmetry constraint: 

Let us first calculate the variational derivatives of A = A(u) with respect to the potentials. 
Taking Gateaux derivative of 4> x = U(u,X)4> at K = (K 1 ,K 2 ) T engenders that 

= (|U + %K 2 + ^'[K})4> + U4/[K\. 
Now by using tp x = —U Tr 4), we obtain 



which leads to 



1 I' 00 dU dU 

*W = -J^. J_J< «*• ^ > * + < «?■ ftT > 



In this way, according to the definition, we get the variational derivatives of A with respect to the 
potentials q, r: 



6X <^ T ,f> 5X <U T ,^> 



** "Too < W r ,ffi > Sr - J_~ < ^ > ^ 

It follows that 



(3.2) 



SX 

~5u E 



1 f°° 
fc J-00 



When the zero boundary conditions: limi^i^-j-oQ (/> — lim| :E |_ i .^ 00 V> = are imposed, we can give out 
a very useful characterization of the variational derivatives of A: 

5X 5X 

tt— = A— , (3.4) 
ou du 

where ^ and 8X/Su are defined by (2.11a) and (3.3), respectively. Introducing N distinct eigenvalues 
Ai, A 2 , • • • , Ajv, we obtain the following spatial and temporal systems 
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1") =V"""(«,A,)f*"Y j = l,2,...,JV, 



IVfe; - "'M^gM-i.V".* (3-56) 



fej =" ^^J' J = l,A---,«, (3 .6a) 



(3.66) 



The compatibility condition of (3.5) and (3.6) is still the mth Kaup-Newell system u tm = K m . For 
(3.5) and (3.6), we may take symmetry constraints 

JG n = J Y^Ej-jr^, or G n = y^^Ej-^j-, Ej = - J (4>ij^ij+r(f>ij^2j -fajipij) dx, n > 0, (3.7) 

where J and G n , n > 0, are defined by (2.11) and (2.10), respectively. In this paper, we are 
interested in the Bargmann constraint (i.e. the first one in (3.7)) 

JG o = J E^« G o=E^- (3.8) 

This kind of constraints requires the G- vector fields to be a potential linear function not including 
any potential differential. The constraint (3.8) yields that 



(S) = (< p 2,gi >,< ap 15 q 2 >) T , 



which gives rise to the following binary nonlinear constraint between the potential and the eigen- 
functions and adjoint eigenfunctions 

q=<AP 1 ,Q 2 >,r=<P 2 ,Q 1 >. (3.9) 

Here for economic writing, we have used the notation: Pj = ((f>n, • • • , 4>iN) T , Qi = (ipn, • • • , ^>iN) T , 
i = 1,2, A = diag(Ai, A2, • • • , Ajv), < y,z >= Ylf=iUj z j f° r two A^-dimensional vectors y = 
(yi,--- ,Vn) T , z = (zi,-- - ,zjv) t G P^. 

We point out that the above explanatory deduction requires Ej / 0, 1 < j < N. But in 
fact, what we need below is just the binary nonlinear constraint (3.9). It itself doesn't touch the 
condition: Ej ^ 0, l<j<N. 

3.2. Nonlinearized Lax pairs and adjoint Lax pairs: 

In what follows, we will denote by P the expression of P(u) under the constraint (3.9). Having 
substituted (3.9) into the Lax pairs and the adjoint Lax pairs: (3.5) and (3.6), we acquire the 
nonlinearized Lax pairs and adjoint Lax pairs 



9ij 




\ 3 r 



q 

-A,- 



-Aj -Aj-r 



r m ~ \m+l 



-q A 

2^i=o °* A 
~ Z^i=o a * A j 



i = i, 2, ■ ■ ■ 
, i = 1,2,- 



•• ,JV; 



771 — I 
j 

m+l — i 



i=o aiA i ~ 2^i=o c i A j 



Em , > 
i=0 6 * A 



m — % 
j 



Em ~ > 
i=0 a iA 



m+1- 




,iV. 



(3.10a) 
(3.106) 
(3.11a) 
(3.116) 



We easily see that the spatial part of the nonlinearized Lax pairs and adjoint Lax pairs, namely, 
the system (3.10) is a finite dimensional system with regard to x, but the temporal parts of the 
nonlinearized Lax pairs and adjoint Lax pairs, namely, the systems (3.11) for n > are all systems 
of evolution equations with regard to t m ,x. Obviously the system (3.10) reads 

( P lx =AP 1 + < AP 1 ,Q 2 > P 2 , 

P 2x =< P 2 ,Qi> AP 1 -AP 2 , 

Qix = -AQi- < P 2 ,Qi > AQ 2 , 

{ Q 2x = - < AP U Q 2 > Qi + AQ 2 , 

and it may be represented as the following Hamiltonian form 



dH _dH_ . 

oQi dPi 



1,2, 



where the Hamiltonian function 



H =< APi,Qi > - < AP 2 ,Q 2 > + < P 2 ,Qi >< AP U Q 2 > 



(3.12a) 



(3.12b) 



In the following, we shall prove that the system (3.10) is an integrable finite dimensional Hamilto- 
nian system in the Liouville sense [17] and that under the control of the system (3.10), the systems 
(3.11) for m > are all transformed into integrable finite dimensional Hamiltonian systems in the 
Liouville sense [17]. 

3.3 An involutive and independent system: 

To show the Liouville integrability of the nonlinearized systems (3.10) and (3.11), we need 
an involutive and independent system to generate integrals of motion for them. We will use the 
nonlinearized Lax pairs and adjoint Lax pairs in the case of zero boundary condition on the eigen- 
functions and adjoint eigenfunctions in order to derive such an involutive and independent system. 
Under the zero boundary condition, the characterization (3.4) asserts that 

' x V„\ / < A m P 2 ,Q 1 > ' 



< A m+1 P 1 ,Q 2 > 



, m > 0, 



(3.13a) 



and further from (2.3), we obtain 

5 = 1, a m = d- l {qZ m -rb m ) = -(< A m P 1 ,Q 1 >-< A m P 2 ,Q 2 >), m > 1. (3.13b) 

It is easy to find that V x = [U, V] still holds. Thus an obvious equality (V 2 ) x = [U, V 2 ] leads to 

Fx = (±trV 2 ) x = ^-(a 2 + \- 1 bZ) = 0. 

This allows us to conclude that F is a generating function of integrals of motion for (3.12). Due to 
that F = J2 n >o P nX~ n , we obtain the following expressions of integrals of motion 

n n — 1 



F n = y^X a i a n-i + ^ C * & ™-1- 



i=0 i=0 

Substitution (3.13) into the above equality gives us the explicit formulas of F n : 



F =1, F x =< AP 1 ,Q 1 > - < AP 2 ,Q 2 > + < P 2 ,Qi >< AP U Q 2 >= H, (3.14a) 

n — 1 1 

Fn=J2j( <Alp i>Qi> ~ < Alp ^ >) (< A n -'Pi, Qi >-< A n ~ l P 2 ,Q 2 >) 



i=l 
n-1 



+ J2< Aip 2,Qi >< A n - l P 1: Q 2 > + < A n P 1 ,Q 1 > -<A n P 2 ,Q 2 >, n>2. 

i=o (3.14b) 

Note that we still have a similar equality V tm = [V^ m \V], m > 0, to (2.12). With the same 
argument, we may find that F = |trV 2 is also a generating function of integrals of motion for 
(3.11). Moreover, a direct calculation can give rise to 

_ d(F m+1 ) _ d(F m+1 ) , . 

Qrt m - Qp — , P*t m - , i - 1,2. (3.15) 



Actually, we have 



(w^ T <^w,^ik^) T '' =1 ' 2 - (316> 



Further it is carried out that 



V dPl ) Z^\ Ci _ a .) f^ j dP 1 \c J -a J 

_ tr W^ h-i\(\ Aj Qi AJ Q2 \ x -a+j) 



i>0 
J>1 



Yl&A'Qi +c i A^Q 2 )\-^\ (3.17) 



i>0 
J>1 



Similarly we can give 



i>l i>0 



j>0 j>l 

It follows from (3.16), (3.17) and (3.18) that 

( 9F m+1 dF m+1 ^T _ (f}(m)\T { Ql 

{ dPi ' dP 2 
A complete similar argument can lead to 



) T = -(V-«™») T (^), m>0. 



The above two equalities enables us to write all the flows of (3.11) in the Hamiltonian form (3.15). 

Now we can easily show the involution F n , n > 1. It is known that for the standard symplectic 
structure on R AN 

to 2 = dQ 1 A dPi + dQ 2 A dP 2 , (3.19) 
the corresponding Poisson bracket is defined as follows 

{f,g} = cu 2 (Idf,Idg), (3.20) 

where Idh expresses a Hamiltonian vector field of a smooth function h on R 4N , defined by Idh Iuj 2 = 
—dh, with J being the left interior product. Because F is also a generating function of integrals of 
motion for (3.11), it is found that 

{F n+1 ,F m+1 } = j^-F n+1 = 0, m,n > 0, (3.21) 

by using the Hamiltonian structure (3.15) of (3.11). This implies that {F n }™ =0 constitutes an 
involutive system with regard to (3.20). It may also directly be shown that {F m ,F n } = 0, m, n > 1. 

Theorem 3.1. Theorem 3.1: Let F k , 1 < A; < N, be defined by 

F k = <t>ikipik + hkip2k, 1 < k < N. 

Then F k , 1 < k < N, F n , n > 1, constitute an involutive system, of which F k , F n , 1 < k,n < N, 
are functionally independent on some region of R 4N . 

Proof: Proof: Noting that we already have {F m ,F n } = 0, m,n > 1, it suffices for the involution 
of Fk, 1 < k < N, F n , n > 1, that we prove {F k ,F n } = 0, 1 < k < N, n > 1. This needs just a 
simple calculation due to the particular form of Fk, 1 < k < N. 
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In what follows, we show that there exists some region of R 4N over which the functions 
Fk, F n , 1 < k,n < N, can become an independent system. Suppose that this result is not true. 
Namely, there doesn't exist any region of R AN over which the functions Fk, F n , 1 < k,n < N, 
could be independent. Hence there exist 2N constants a k , (3 n , 1 < k,n < N, satisfying 



N 



so that we have for all points (Pi , P 2 , Qj , Q 2 ) m R AN 



N 

k=i 



T 



(3.22) 



From the above equality, we would like to derive that all constants a k , (3 n , 1 < k,n < N, must be 
zero. We will utilize the following equalities: 



dF k 



dPi 
dF 1 



(0,--- ,0,V, fc ,0,--- ,0) T , i = l,2, l<k<N, 



Pi=P 2 =0 



k-i 



dP x 

dF n 



dF t 



Pi=P 2 =0 



dP 1 



= A n Qi 



dP 2 
dF 



N-k 



Pi=P 2 =0 



-AQ 2 , 



Pi=P 2 =0 



dP 2 



= -A n Q 2 , n>2. 



Pi=P 2 =0 



The equality (3.22) upon choosing Pi = P 2 = engenders 



= Y J PnX n k =-Y,^K, l<k<N. 



n=l 



n = l 



Therefore we have J2n=i Pn^k = 0, 1 < k < N. Let /(A) = Yln=i Pn^ n , which is a polynomial 
of degree N. Since it possesses N distinct roots Ai, A2, • ■ ■ , Xn, we have /(A) = 0. This means 
that f3 n = 0, 1 < n < N, and further we obtain a fe =0, l^A;^^. In this way, we arrive at 
all zero constants, which contradicts to our assumption. Hence there must exists at least a region 
f2 C R 4N over which the 2N 1-forms dF k: dF n , 1 < k, n < N, are linearly independent. The proof 
is finished. 

The theorem provides us with an involutive and independent system of polynomials, which is 
completely new as far as we know. 

3.4 Liouville integrability of the nonlinearized Lax pairs and adjoint Lax pairs: 

We now proceed to show the Liouville integrability of the nonlinearized Lax pairs and adjoint 
Lax pairs, i.e. (3.10) and (3.11). We point out that the above manipulation with the explicit 
expression (3.15) of (3.11) requires the zero boundary conditions: lim| a ,|_ >+00 Qi = \im\ x \^ +ryo Pi = 
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0, i = 1,2. Now let us deal with a general case, i.e. the case where the zero boundary conditions 
are not imposed. In this case, it follows from 4> x = U(u,X)(f>, if) x = —U(u, \) T ip that 

where / is an integral of motion of <fi x = U(u, X)(j), ip x = —U(u, X) T tp. By applying (3.23) n times, 
we easily work out that 

(ty^HA^x^^ <-> 

where /o = 1 an d 1%, 1 < i < m, are all integrals of motion of (3.10). Further we can compute that 
a m = <9 -1 (gCm -rb m ) 

1 m — l 

=- h(< A m -*P U Qi > - < A m ~ i P 2 , Q 2 >) + T m , m > 1, (3.25) 

i=0 

where T m is an integral of motion of (3.10), too. The latter two equalities in (2.3) lead to that 
T m = —I m , m > 1. Therefore by using (2.5), we obtain 



m— 1 



h{< A m ~ l P u Q 1 >-< A m ~ l P 2 ,Q 2 >) - 2I m 

i=0 

m—l i—1 

= - E 4 [E J ^< ^ _fcp i^i > - < At ~ kp ^ >) - 2/ * 

i=l k=0 
m — i — 1 

[J2 h{< A m ~ l ~ l P 1 , Qi> - < A m ~ l ~ l P 2 ,Q 2 >) - 2/ m . 

1=0 

m—l i m — l — i 

~ E E /fc < ^ +1_fc ^i^2 > E /; < ^ m " 1_l "^2,Qi >, m > 2. (3.26) 

i=l fc=0 2=0 

By interchanging the summing in the above equality and noting the expressions of F^s, a direct 
but quite lengthy calculation may give rise to 

1 1 m_1 

Im = 2 E IkhF m -( k +i) - 2 E ^rn-i, m > 2. (3.27) 

k+l<m — 1 i=l 
fc,Z>0 

This equality and Ji = determined by cti = — | gr, readily yield the following expressions of 



/m = Ed» E *ii ™>1- (3-28) 



n=l ii H \-i n =m 
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Here the constants d n are denned recursively by 

n — 1 n — 1 



13 1 1 C 1 ' ^ 

di = -, d 2 = g, = d n -i + 2 d s dn-s-i - 2 E d sdn-s, n > 3. (3.29) 

s=l s=l 

At this point, we can give the expressions of a m , 6 m , c m in terms of F„ and further a direct 
computation may show that the temporal parts (3.11) of the nonlinearized Lax pairs and adjoint 
Lax pairs are represented as the following Hamiltonian systems 

^=-^,^ = ^,-1.2 (3.30.) 
with the Hamiltonian functions 

F - = E^fl E F il ---F i „ +1 (do = l), m>0. (3.30b) 

n=0 iiH \-i n + 1 =m+l 

,in + l>l 

The Hamiltonian phase flows g^ m of the systems (3.30) commute with each other due to the 
commut ability of F n . 

Because H, H m , m > 1, are all functions of F n , n > 1, we can now conclude the Liouville 
integrability on the nonlinearized Lax pairs and adjoint Lax pairs. 

Theorem 3.2. Theorem 3.2: The spatial part (3.10) of the nonlinearized Lax pairs and adjoint 
Lax pairs and the temporal parts (3.11) of the nonlinearized Lax pairs and adjoint Lax pairs under 
the control of the spatial part are all finite dimensional Liouville integrable Hamiltonian systems 
with the corresponding Hamiltonian functions H and H m defined by (3.13b) and (3.30b). Moreover 
they all possess involutive and independent integrals of motion: F^, F n , 1 < k,n < N. 

3.5. Involutive solutions to Kaup-Newell hierarchy: 

The above manipulation also allows us to establish the following involutive solutions of the 
mth Kaup-Newell system u tm = K m : 

( q(x,t m ) =< A-Sii9}i m P\ (0,0), SH9n m Q2 (0, 0) > (3.31a) 
1 r(x,t m ) =< ff^ m A(0,0), 5 ^ m Q 2 (0,0) > (3.316) 

where g^, g l ff m are the Hamiltonian phase flows [17] associated with the Hamiltonian systems 
(3.12) and (3.30), respectively, but Pi(0, 0), Qi(0,0), i = 1,2, may be arbitrary. This provides us 
with a way to present solutions to Kaup-Newell systems by solving two separate finite dimensional 
integrable systems with the variables x and t m . In fact, in view of the Liouville integrability of the 
nonlinearized Lax pairs and adjoint Lax pairs, this kind of involutive solutions to integrable systems 
not only gives rise to the interrelation between 1 + 1 dimensional integrable systems and finite 
dimensional integrable systems, but also exhibits the characteristic of integrability by quadratures 
for integrable systems in 1 + 1 dimensions. 
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4. Conclusions and remarks 



We have introduced a different spectral problem from the original one for Kaup-Newell systems 
and discussed the corresponding Bargmann symmetry constraint. We have also exhibited a new ex- 
plicit involutive system {Fk, 1 < k < N, F n , n > 1} on the symplectic manifold (R AN ,u> 2 ), of which 
{Ffc, F n , 1 < k, n < N} constitute an independent system. The Lax pairs and adjoint Lax pairs 
are nonlinearized into a hierarchy of finite dimensional Liouville integrable Hamiltonian systems 
and Bargmann symmetry constraint leads to a sort of involutive solutions to Kaup-Newell systems, 
which may be presented by solving two commutative finite dimensional Hamiltonian systems, i.e. 
the spatial part and the temporal part of nonlinearized Lax pairs and adjoint Lax pairs. 

We remark that under the reduction Qi = —P2, Q 2 = Pi, the involutive system {Fnj^i 
defined by (3.14) are reduced to 

F l = - 2 < AP 1 ,P 2 >-< AP 1 ,P 1 >< P 2 ,P 2 >, 

n 

F n =J2{< Aip ^ p 2 >< A n ~ i P 1 ,P 2 >-< A t P 1 ,P 1 >< A n ~ i P 2 ,P 2 >) 

i=0 

- 2 < A n P 1 ,P 2 >-< P U P 2 >< A n P u P 2 > + < P 1 ,P 1 >< A n P 2 ,P 2 >, n > 2. 

This involutive system corresponds to mono-nonlinearization for the Kaup-Newell systems (2.10) 
(see Ref. [18]). It involves 2N dependent variables but not AN dependent variables. 

It should be noted that binary nonlinearization procedure is theoretically very systematic (refer 
to [19] for more information). It paves a way to show Liouville integrability for both the spatial 
part (sometimes called constrained flows or restricted flows) and the temporal parts of nonlinearized 
Lax pairs and adjoint Lax pairs, which is short of discussion in the works of restricted flows (see 
for instance [8]). We only need the distinction of Ai, A2, • • • , Xn for proving Liouville integrability 
but may have the condition Ylf =1 Aj = 0, which is not allowed in the case of restricted flows. The 
resulting binary nonlinear constraints involve two classes of eigenfunctions and thus might greatly 
extend explicit finite dimensional integrable Hamiltonian systems in the Liouville sense. 
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